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We present an observation for ZZ — > £'^£~£ ~ {£, £ = e or n) production in pp collisions at a 
center-of-mass energy of = 1-96 TeV. Using 1.7 fl)~^ of data collected by the DO experiment at 
the Fermilab Tevatron Collider, we observe three candidate events with an expected background of 
0.141q q2 events. The significance of this observation is 5.3 standard deviations. The combination 
of DO results in this channel, as well as in ZZ —^ £'^£~vi>, yields a significance of 5.7 standard 
deviations and a combined cross section of a{ZZ) = 1.60 ± 0.63 (stat.)^^ !^ (syst.) pb. 
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Studies of the pair production of electroweak gauge 
bosons provide an interesting test of the electroweak 
theory predictions [J]. In contrast with other diboson 
processes, Z boson pair production {ZZ) does not in- 
volve trilinear gauge boson interactions within the stan- 
dard model (SM). The observation of an unexpectedly 
high cross section could indicate the presence of anoma- 
lous ZZZ or ZZ^ couplings The SM prediction for 
the total ZZ production cross section in collisions 
at the Fermilab Tevatron Collider at %fs — 1.96 TeV is 
a{ZZ) = 1.4 ±0.1 pb The requirement of leptonic Z 
boson decays reduces the observable cross section, mak- 
ing its measurement rather challenging. The accumu- 
lation of integrated luminosities in excess of 3 fb~^ at 
the Fermilab Tevatron Collider and the development of 
highly optimized event selection criteria has now made 
possible the direct observation of ZZ production. 

Previous investigations of ZZ production have been 
performed both at the Fermilab Tevatron and the 
CERN e+e" (LEP) d CoUiders. The DO collaboration 
reported a search for ZZ t^l^ I ^ {£, £ — e or /i) 
with 1 fb~^ of data that provided a 95% C.L. limit of 
a(ZZ) < 4.4 pb The CDF collaboration reported 
a signal for ZZ production with a significance of 4.4 
standard deviations from combined ZZ £^£~£ ^£ ~ 
and ZZ — > £'^ £~vv searches, and measured a production 
cross section of a[ZZ) = 1.4^°;^ pb @. 

In this Letter, we present a search for Z boson pairs 
where the Z bosons have decayed to either electron or 
muon pairs, resulting in final states consisting of four 
electrons (4e), four muons (4/i) or two muons and two 
electrons (2/i2e). Data used in this analysis were col- 
lected with the DO detector at the Fermilab Tevatron 
ColUder at Vs = 1-96 TeV between June 2006 and May 
2008. The integrated luminosities 01 for the three an- 
alyzed channels are about 1.7 fb^-'^. This result is later 
combined with that from an earlier similar analysis Q us- 
ing data collected from October 2002 to February 2006 
and corresponding to an integrated luminosity of 1 fb~^. 

The DO detector Q consists of a central tracking sys- 
tem, comprised of a silicon microstrip tracker (SMT), 
and a central fiber tracker (CFT), providing coverage to 
pseudorapidity [yyl < 3 both located within a 2 T 
superconducting solenoidal magnet. Three liquid argon 
and uranium calorimeters provide coverage to \r\\ < 4. 
Electromagnetic objects are well reconstructed in the re- 
gions of the central calorimeter (CC) with coverage to 
I77I < 1.1 and the end calorimeters (EC) with coverage of 
1.5 < I77I < 3.2. A muon system surrounds the calorime- 
try, consisting of three layers of scintillators and drift 
tubes and 1.8 T iron toroids, with a coverage of \r]\ < 2. 

This analysis employs a combination of single and di- 
electron triggers for the 4e channel. Similarly, single and 
dimuon triggers are used for the Afj, channel. The 2/z2e 
channel uses a combination of all these triggers, and ad- 
ditional specific electron-muon triggers. The triggering 



efficiency for events with four leptons having high trans- 
verse momentum (pt) that satisfy all offline selection re- 
quirements exceeds 99%. 

For the 4e channel, we require four electrons with or- 
dered transverse energies Et > 30, 25, 15, and 15 GeV, 
respectively. Electrons must be reconstructed either in 
the CC region or in the EC region, and be isolated from 
other energy clusters in the calorimeter. Electrons in the 
CC are required to satisfy identification criteria based 
on a multivariate discriminant derived from calorimeter 
shower shape and a matched track reconstructed in the 
SMT and CFT. Electrons in the EC are not required to 
have a matched track, but must satisfy more stringent 
shower shape requirements. At least two electrons must 
be in the CC region. With no requirement applied on 
the charge of the electrons at this stage to increase se- 
lection efficiency, three possible ZZ combinations can be 
formed for each 4e event. Events are required to have 
a solution for which one ee combination has an invari- 
ant mass > 70 GeV and the other > 50 GeV. Finally, 
events are split into three categories, depending on the 
number of electrons in the CC region. Subsamples with 
two electrons, with three electrons and with four or more 
electrons in the CC are denoted as 4e2C, 4e3c, and 4e4c, 
respectively. The three exclusive subsamples contain sig- 
nificantly different levels of background contamination 
and thus the separation of the subsamples provides more 
sensitivity to the search. 

For the 4/i channel, each muon must satisfy quality cri- 
teria based on scintillator and wire chamber information 
from the muon system, and have a matched track in the 
central tracker. We require that the four most energetic 
muons have ordered transverse momenta pT > 30, 25, 
15, and 15 GeV, respectively. At least three muons in 
the event must be isolated, each passing a requirement 
of less than 2.5 GeV of transverse energy deposited in 
the calorimeter in the annulus 0.1 < AR < 0.4 centered 
around the muon track [13]. Finally, the muon IS re- 
quired to be well reconstructed and to originate from the 
primary event vertex. Of the three possible ZZ combina- 
tions per event that can be formed without considering 
muon charge at this stage, a solution is required where 
one fifi combination has an invariant mass > 70 GeV and 
the other > 50 GeV. 

For the 2/i2e channel, the two most energetic electrons 
and muons in an event must have Et{pt) > 25, 15 GeV. 
All muons and electrons must satisfy the single lepton 
selection criteria defined for the 4e and 4/i final states. 
In addition, electrons and muons are required to be spa- 
tially separated by AR > 0.2 to remove Z fifj, back- 
ground with muons radiating photons giving events with 
two muon and two trackless electron candidates. At least 
one muon must satisfy the same calorimeter isolation re- 
quirement imposed in the 4/i final state. A solution is 
required where one pair of same flavor leptons has an 
invariant mass > 70 GeV, and the other > 50 GeV. Fi- 
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nally, events are split into three categories depending on 
the number of electrons in the CC region. Subsamples 
with no electron, with one electron and with two or more 
electrons in the CC are denoted as 2/i2eocj 2/i2eic a-nd 
2/i2e2Cj respectively. As in the 4e channel, these sub- 
samples have significantly difi^erent levels of background 
contamination. 

A Monte Carlo (MC) simulation is used to determine 
the expected number of signal events in each subchan- 
nel. The small contribution from ZZ events with at 
least one Z boson decaying into tau pairs is also in- 
cluded in the signal. Simulated events are generated us- 
ing PYTHIA [ll| and passed through a detailed GEANT- 
based simulation IJ] of the detector response. Differ- 
ences between MC simulation and data in the recon- 
struction and identification efficiencies for electrons and 
muons are corrected using efficiencies derived from large 
data samples of inclusive Z II {I — e, or [i) events. 
The systematic uncertainty in the signal is dominated by 
the uncertainty in the theoretical cross section (6.25%), 
the uncertainty on the lepton identification and recon- 
struction efficiencies (« 4% for the 4e and 4/i subchannels 
and « 2.5% for the 2/Lt2e subchannels) and a 6.1% un- 
certainty on the luminosity measurement Q ■ Additional 
smaller sources of systematic uncertainty arise from en- 
ergy and momentum resolutions and MC modeling of the 
signal process. 

Backgrounds to the ZZ signal originate from top quark 
pair (if) production and from events with W and/or Z 
bosons that decay to leptons and additional jets or pho- 
tons. The jets can then be misidentified as leptons or 
contain true electrons or muons from in-flight decays of 
pions, kaons, or heavy-flavored hadrons. 

The background from tt production is estimated from 
MC, assuming the cross section of a{t£) = 7.9 pb 
for a top quark mass of 170 GeV. The systematic un- 
certainty includes a 10% uncertainty on a{ti), as well as 
contributions from the variation in cross section and ac- 
ceptance originating from uncertainties on the mass of 
the top quark. 

To estimate the misidentified lepton background, we 
first measure the probability for a jet to produce an elec- 
tron or muon that satisfy the identification criteria from 
data using a "tag and probe" method The probabil- 
ity for a jet to mimic an electron, parameterized in terms 
of jet Et and ry, is equal to 4 x 10"'' for the case of CC 
electrons with a matched track and 5 x 10~^ in the case 
of EC electrons for which no track matching is applied. 
The probability for a 15 GeV (100 GeV) jet to produce 
a muon of > 15 GeV is 10~''(10~^) without requiring 
muon isolation, and it is 10^^(10^'') when the muon is 
required to be isolated. A systematic uncertainty of 30% 
on the jet-to-lepton misidentification probabilities is es- 
timated by varying the selection criteria of the control 
samples used. 

The probabilities for jets to be misidentified as elec- 
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FIG. 1: Distribution of four lepton invariant mass in data, 
expected signal, and expected background. 



trons are then applied to jets in eee-|-jets and fifie+jets 
data to determine the background to the 4e and 2/i2e 
channels, respectively. This method takes into ac- 
count contributions from Z+jets, Z+j+iets, WZ+jets, 
WW+jets, W^-t-jets, and events with > 4 jets. However, 
this method double counts the contribution from Z+jets. 
A correction is measured, amounting to « 20%, to cor- 
rect for the double counting. The probabilities for jets 
to contain a muon are applied to jets in /i/^+jcts data 
to determine a background estimate for the 4/i channel. 
Systematic uncertainties on this background arise from 
the 30% uncertainty in measured misidentification rates, 
and from the limited statistics of the data remaining in 
the samples after selection. 

Table U summarizes the expected signal and back- 
ground contributions to each subchannel, as well as the 
number of candidate events in data. The total signal 
and background expectations are 1.89 ± 0.08 events and 
0.14to;o2 events, respectively. We observe a total of three 
candidate events in the data, two in the 4e4c subchannel 
and one in the 4/i subchannel. Table [III summarizes some 
of their kinematic characteristics. The quoted dilepton 
masses in the table correspond to the one out of the three 
possible combinations having opposite charge within the 
pairs and having a dilepton mass closest to that of the Z 
boson. Figure [T] shows the distribution of the four lepton 
invariant mass for data and for the expected signal and 
background. 

We extract the significance of the observed event dis- 
tributions using a negative log-likelihood ratio (LLR) 
test-statistic [ij]. As input, we use the expected yields 
(number of events) from signal and background, sepa- 
rated into the seven subchannels compared to the ob- 
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TABLE I: The integrated luminosity, expected number of signal {Z/'y* Z/'y") and background events [tt and Z(7)+jets which 
includes all W / Z /'y+]ets contributions), and the number of observed candidates in the seven ZZ I'^l'l ~ subchannels. 
Uncertainties reflect statistical and systematic contributions added in quadrature. 



Subchannel 


4e2C 


4e3C 


4e4c 


4m 


2^2eoc 


2^i2eic 


2^i2e2C 


Luminosity (fb~^) 


1.75 ±0.11 


1.75 ±0.11 


1.75 ±0.11 


1.68 ±0.10 


1.68 ±0.10 


1.68 ±0.10 


1.68 ±0.10 


Signal 


0.084 ± 0.008 


0.173 ±0.015 


0.140 ±0.012 


0.534 ± 0.043 




0.352 ± 0.040 


n trtro + 0.045 


2'(7)+jets 


U.UOU_f, QQg 


o.oi8tE::r 




0.0003 ± 0.0001 


03+° "^ 


0.05 ±0.01 


0.008„(, 


ti 










0012+" °°^'^ 

U.UUiZ_Q QQQ9 


0.005 ± 0.002 


0007+°-'"'°^ 

U.UUU/ _o 0005 


Observed events 








2 


1 












TABLE II: Characteristics of the observed candidate events. 
r\ and values are measured relative to the location of the pp 
collision. Mu is the mass of the lepton pair. 



4e 

candidate 1 




4 


e+ 

^2 


63 


^4 


PT (GeV) 


107 


59 


52 


16 


■n 


0.66 


0.25 


-0.64 


-0.85 




4.10 


1.08 


0.46 


2.62 


Mu (GeV) 


6^64 
89 ±3 


61 ±2 


4e 

candidate 2 




4 


4 


63 


64 


PT (GeV) 


83 


75 


35 


26 




0.64 


0.40 


0.85 


1.17 




6.16 


3.80 


3.83 


1.40 


Mu (GeV) 


99 ±3 


6^64 
90 ±4 


4m 
candidate 




4 


M2 


M3 




PT (GeV) 


115 


77 


42 


24 




0.04 


-1.01 


0.77 


-1.93 




1.69 


4.26 


5.29 


0.36 


Mu (GeV) 




M2 mI 

90t^2 



served yields. The modified frequentist method returns 
the probability (p-value) of the background-only fluctu- 
ating to give the observed yields or higher. In 5x10^ 
background pseudo-experiments, we find 213 trials with 
an LLR value smaller or equal to that observed in data. 
This gives a p- value of 4.3x10^^ which corresponds to a 
5.3 standard deviation (cr) observed significance (3.7cr ex- 
pected). The probability for the signal plus background 
hypothesis to give less signal-like observations than the 
observed one is 0.87. A correction factor of 0.93, de- 
rived using PYTHIA, is used to convert the measured cross 
section for {Z/j*){Z/^*) into that for ZZ production. 
Minimizing the LLR function we obtain a cross section 
of a{ZZ) = 1.75+J;86 (stat.) ± 0.13 (syst.) pb for this 
analysis. 

This result is combined with the results from an in- 
dependent ZZ I'^t^vv search [l^, and the previ- 
ous (Z/7*)(Z/7*) ^ analysis which used 
a separate data sample with a looser mass requirement 



M{U) > 30 GeV. The earlier search contributes no signal 
events and we have scaled its background estimate to the 
tighter kinematic range used in the recent analysis. The 
combination of the three analyses is performed taking 
into account the correlations of systematic uncertainties 
between subchannels and among analyses. The resulting 
p-value is 6.2x10^^, and the significance for observation 
of ZZ production increases to 5.7(7 (4.8(7 expected). The 
probability for the signal plus background hypothesis to 
give less signal-like observations than the observed one 
is 0.71. We therefore report the observation of a ZZ sig- 
nal at a hadron collider with a combined cross section 
of (7(ZZ) = 1.60 ±0.63 (stat.)+° :[7 (syst.) pb, consistent 
with the standard model expectation. 
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